Kinematics of geodesic flows in stringy black hole backgrounds 



Anirvan Dasgupta* 
Department of Mechanical Engineering and Centre for Theoretical Studies 
Indian Institute of Technology, Kharagpur 721 302, India 

Hemwati Nandan''' 
0> ' Centre for Theoretical Studies 

o ■ 

O ■ Indian Institute of Technology, 

CN 

C5 \ Kharagpur 721 302, India 

CSI ' Sayan Kar^ 

Department of Physics and Centre for Theoretical Studies 

qh| Indian Institute of Technology, Kharagpur 721 302, India 

in . 

m 
> 

o 

cn 

o 

oo 

o 

> 

X 
J3 



* Electronic address: amr@mech.iitkgp.ernet.in 
^ Electronic address: hnandan@cts.iitkgp.ernet.in 
^ Electronic address: sayan@cts.iitkgp.ernet.in 

1 



Abstract 

We study the kinematics of timelike geodesic congruences in two and four dimensions in space- 
time geometries representing stringy black holes. The Raychaudhuri equations for the kinematical 
quantities (namely, expansion, shear and rotation) characterising such geodesic flows arc written 
down and subsequently solved analytically (in two dimensions) and numerically (in four dimensions) 
for specific geodesies flows. We compare between geodesic flows in dual (electric and magnetic) 
stringy black hole backgrounds in four dimensions, by showing the differences that arise in the 
corresponding evolutions of the kinematic variables. The crucial role of initial conditions and the 
spacetime curvature on the evolution of the kinematical variables is illustrated. Some novel gen- 
eral conclusions on caustic formation and geodesic focusing are obtained from the analytical and 
numerical findings. We also propose a new quantifier in terms of the time (affine parameter) of 
approach to a singularity, which may be used to distinguish between flows in different geometries. 
In summary, our quantitative findings bring out hitherto unknown features of the kinematics of 
geodesic flows, which, otherwise, would have remained overlooked, if we confined ourselves to only 
a qualitative analysis. 

PACS numbers: 04.20.Cv, 83.10.Bb, 04.40.-b, 97.60.Lf 
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I. INTRODUCTION 



The kinematics of Efeodesic congruences is characterised by three kinematical quantities: 
isotropic expansion, shear and rotation (henceforth referred as ESR) QflSflH,!^. The 



evolution of these quantities along the geodesic flow, are obtained from the Raychaudhuri 
equations These equations are derived by relating the evolution of the 

deformation (or deviation) vector between two neighbouring geodesies (expressed in terms 
of the ESR variables) to the curvature of the space/spacetime. This is already a well-studied 
subject (see Q and the references therein). 

In its original incarnation, the Raychaudhuri equation provided the basis for the descrip- 
tion and analysis of spacetime singularities in gravitation and cosmology Q]. For example, 
the equation for the expansion and resulting theorem on geodesic focusing is a crucial ingre- 
dient in the proofs of Penrose-Hawking singularity theorems js, 9|. However, these equations 
have a much wider scope in studying geodesic as well as non-geodesic flows in nature which 
may possibly arise in diverse contexts (see for some open issues). We have recently used 
these equations to investigate the kinematics of flows on flat and curved deformable media 

nn 

(including elastic and viscoelastic media) in detail |10l . . 

In this article, we attempt to understand the kinematics of geodesic flows in the presence 
of spacetime geometries representing black holes. Though, it is true that the Raychaudhuri 
equations have been around now for more than half a century, we are not aware of any 
attempt at a complete study of its solutions and the dependence of the ESR variables on the 
initial conditions imposed on them (however, see for a recent work). A slight subtlety 
may be noted here. The Raychaudhuri equations are for the ESR variables, but they also 
involve the tangent vector fleld (denoted as u\ later). Thus, unless one knows the solutions 
for the (i.e. the flrst integrals of the geodesic equations) one cannot proceed towards 
solving for the ESR. Alternatively, one can try to flnd solutions for the full set of variables, 
i.e. as well as the ESR, by imposing initial conditions on all of them and evolving the full 
system of equations along the flow. This approach enables us to obtain the tangent vector 
fleld as well as the ESR simultaneously. In our work, we adopt this method primarily for the 
four dimensional cases where we are unable to solve for the geodesies or the analytically. 

As an aside, it may be noted, that the study of accretion of matter, or evolution of space- 
time deformations near a black hole has been a fascinating topic of study in classical general 
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relativity. Such processes can perhaps be investigated using the Raychaudhuri equations. 
The study of kinematics of non-spacehke congruences may further provide useful insights on 
the observable features in a given spacetime. We take this as a background motivation for 
our study. 

The non-trivial geometry of different types of black holes is of crucial interest in general 
relativity and for our purpose, we will consider some typical two and four dimensional 
stringy black hole spacetimes. The two dimensional black hole geometry we work with was 
first obtained in the context of string theory j3|-jl6| in the nineties by Mandal et. al. fl?! . 
Exact solutions for geodesies and geodesic deviation in this two dimensional stringy black 
hole background have already been studied earlier In four dimensions, the simplest 

eternal black hole geometry is, of course, the Schwarzschild. Rather than working with just 
the Schwarzschild alone, we consider the variations of the Schwarzschild which have arisen 

m 



in the context of string theory 
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2l| . All the four dimensional geometries we choose 



to work with have a Schwarzschild limit (obtainable by setting a parameter in the line 
element to zero). Further, the validity of the various (weak, strong, null averaged) energy 
conditions for the matter that threads such stringy black hole spacetimes have also been 



investigated in detail in order to understand the nature of geodesic focusing 



22|. 



Our article is organised as follows. We first review the background spacetimes (Section 
II) and then (Section III) derive the evolution (geodesic and Raychaudhuri) equations for 
the two dimensional stringy black hole metric. In this 2D case, the Raychaudhuri equa- 
tion is just the equation for the expansion scalar. We obtain analytical solutions for the 
expansion scalar, based on which the role of initial conditions on caustic formation and 
geodesic focusing/defocusing is investigated and analysed. Subsequently (Section IV), we 
turn to actual four dimensional solutions in dilaton-Maxwell gravity. We analyse the nature 
of deformations for the geodesic flows in the Garfinkle-Horowitz-Strominger (GHS) electric 
and magnetic (dual) solutions [l9] and compare our results with those for Schwarzschild 
geometry (obtainable from the GHS metrics by setting a stringy parameter to zero). Here, 
the kinematics of deformations is studied on the two dimensional equatorial plane. In the 
absence of analytical solutions, we solve the geodesic and Raychaudhuri equations numeri- 
cally under differing initial conditions on the associated variables. The generic features of 
the evolution of deformations are brought out. The influence of the gravitational field on the 
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evolution of the ESR is discussed and the effect of curvature is understood. Finally (Section 
V), we conclude our results and suggest some relevant future work. 



II. THE STRINGY BLACK HOLE SPACETIMES 



In this section, we quickly recall a well-known solution obtained in the context of two 



dimensional low energy effective string theory 



2l| . The line element which we mention 



below solves the so-called /3-function equations for the string cr-model 2l|]. It is known that 



with appropriate methods of compactification one can obtain effective equations in two as 



well as other higher dimensions. These equations, in the simplest scenario, are for t 



le metric 



field and the dilaton field. We quote below the line element in two dimensions [13|, Il7l |. 



where, m and k, (with dimensions of length square) are the mass and central charge param- 
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eters, respectively, and are linked with the concepts in two dimensional string theory 
which we do not bother about here. 

In 3 + 1 dimensions, we also have asymptotically flat solutions representing black holes in 
dilaton-Maxwell gravity. Such solutions, due to Garfinkle, Horowitz and Strominger jioj l. 
represent electric and dual magnetic black holes [2^ . The spacetime geometry of these line 
elements are causally similar to Schwarzschild geometry. The metric for the black hole with 
electric charge is given as, 

ds' = dt^ + + dnl (2.2) 

m sinh a \ I -L „ J 



1 + 

where dQ^ = {dip"^ + sin^ ip dcj?) is the metric on a two dimensional unit sphere and a is 
a parameter related to the electric charge. Further, the dual (magnetic) metric of (12. 2p is 
given as follows ^|, 

(1-81) (1 _ 22) (1 _ Qi) 

where Q is the magnetic charge of the black hole. In the respective limit of a = or Q = 0, 
we have Schwarzschild geometry in both the cases. We shall investigate the ESR variables 
for geodesic flows in each of the above two four dimensional metrics and compare our results 
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for the electric and magnetic solutions with those for the Schwarzschild. 

For a general and compact representation of the above spacetimes, one can use a generic 

line element in 3 + 1 dimensions as follows, 

ds^ = -X{r) + Y{r) dr^ + dfll (2.4) 

for specific choices of X{r) and Y{r). The general structure of the geodesic equations for 
the line element (12. 4p are given by, 

t+^rt = 0, (2.5) 
/X'(r)t2 + y'(r)f2 _ 2^ J,2 _2r sin^TA 02 \ ^ ^ 

' + — ) = «■ 

■■ 2 • 

'?/' H — rip — cos Tp sin t/; 0^ = 0, (2-7) 
r 

.. 2 • 

0+-r0 + 2cot^/'^/'0 = O, (2.8) 
r 

where the prime denotes the differentiation with respect to r. Using the functional forms 
of X{r) and Y{r) (corresponding to the line elements given by (12.21) and (12. 3p . in equations 
(12. 5p and (12. 6p . one can easily obtain the set of geodesic equations for the particular cases of 
the electric and magnetic stringy black holes. For the line element (12. ip in 1 + 1 dimensions, 
the term dVL\ is absent in (12. 4p and the geodesic equations are given by (12.51) and (12. 6p 
(without the terms if) and 0). 

III. KINEMATICS OF DEFORMATIONS IN 1 + 1 DIMENSIONS 
A. Kinematic variables 

The evolution of space-like deformations in a two dimensional geodesic congruence is 
captured through the evolution of the geodesic deviation vector (where i = 1,2) on a 
space-like hypersurface. These deformations can be described in terms of a second rank 
tensor Bj = V jU^ [11, llOj, which governs the dynamics of the congruence. The second order 



derivative of the vector with respect to an affine parameter is given as follows 10|, lUl , 

C = {Bi + B\B'^e . (3.1) 
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The evolution tensor Bj in equation (13. ip is usually decomposed into irreducible parts sig- 
nifying the expansion (scalar 6), shear (trace-free tensor a*) and rotation (antisymmetric 
tensor Uj). Since MjCr] = (space-like deformations), and using the zero trace property, i.e., 
a\ = 0, leads to ctq = cr\ = ctq = 0. Similarly, the rotation tensor uj^j also satisfies Mjcu^ = 
which leads to uj\ = = 0. Therefore, the evolution tensor can be expressed only in terms 
of the expansion scalar, in the following form, 

B'j = eh), (3.2) 

where the projection metric is defined as h^j = Uj. Here, Ui is a time-like vector field. 

B. The evolution equations 

The evolution equations for a congruence of time-like geodesies for the present case consist 
of the Raychaudhuri equation for the expansion scalar and the geodesic equations derived 
for a particular metric. In order to derive the Raychaudhuri equation for the expansion 
scalar, we first write down the second derivative of the deformation vector in the following 
form, 

c = -R\,yu^e- (3.3) 

1. Raychaudhuri equation for expansion scalar 

Using the equation (13. 2p and (13.30 in the equation (13. ip . one can now obtain the Ray- 
chaudhuri equation for the expansion scalar as given below, 

e + e"^ = -R'n^u^vJ^ = -R imu^u"". (3.4) 

The general form of equation (13. 4p for the metric (12.11) can also be written as follows: 

e + e'-^ = o, (3.5) 

where the Ricci scalar R = Am/ nr. It may be noted that for r —>■ oo, the equation (13. 5p 
reduces to that in flat space without shear and vorticity. This evolution equation (13. 5p 
along with the geodesic equations corresponding to the metric (12. ip . form a complete set of 
equations required for studying the kinematics of deformations of geodesic congruences. 
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2. First integrals of geodesic equations 



The first integrals of tlie geodesic equations in 1 + 1 dimensions can be obtained as, 



E 

2(1 



i = ... m, ^ (3.6) 



r 



= - \{E'^ -4:)r'^ + 4:mr] = -V2ir), (3.7) 

K 

wliere V2(r) is an effective potential. It may be noted that the constraint Qapu'^u^ = — 1 
(time-like geodesies) is used to obtain (13. 7p . The different choices for the constant of motion 
E result in the different behaviour of the effective potential. One can have a harmonic, or an 
inverted harmonic oscillator corresponding to E"^ < 4, or E"^ > 4, respectively, while E"^ = 4 



results in a linear potential with negative slope 



18j. We will now solve the Raychaudhuri 



equation for the expansion scalar for all the choices of E mentioned above. 



C. Exact solution for expansion scalar 

The equation (13.51) can be solved for the above-mentioned three different cases by inte- 
grating equation (13. 7p once (see [ISl]), and then using r in equation (13.51) . 
Case (A): E^ <A. 



The equation (13. 5p for the expansion scalar with E"^ < A reads: 

^(A) + ^2^A) -2sec2A = 0, (3.8) 

where we have used the scaling {A,^} = [(4 — £'^)/4k]^/^{A, 6'}. The solution of equation 
(13. Sp is then given by, 

- ^ = tanA + (D, + A)sec^A 

^ ' l + (Di + A)tanA ^ ' 

where Di is an integration constant which can be given in terms of the initial conditions as 
follows, 

^ (1 - Apgo) tan Aq + Ap sec^ Ap - Op 
Oq tan Ap — sec^ Ap 

Case (B) : E^ = A. 

The equation for expansion scalar with i?^ = 4 reads: 

e{x) + e\x)-^ = o, (3.11) 



A 
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and the solution of equation (13.111) can be given as follows, 



where the integration constant D2 is given as, 

(1 + Apgo) 

- -A3(Vo-2)- ^'-''^ 

Case (C) : E^>A. 

The equation for the expansion scalar for this case reads: 

^(A) + ^2(A)-2cosech^A = 0, (3.14) 

where {A,^} = [(E^ - 4)/4k]1/2{A, 0}. The solution of equation fIXTD is, 

- - coth A - (D3 + A) cosech^A 

= (D3 + A)cothA-l • ^'-''^ 



^3 = ^ Y , T 2Y • (3-16) 



where the integration constant D-^ is given as, 

Aq^o) coth Aq — Aq cosech^Ao + Oq 
9q coth Aq + cosech^Ao 

It is in order to mention here that, following an altogether different approach, one can 
calculate the expansion scalar ^ as a function of r from the expressions of the first integrals 

= {i, r) (the velocity field) given by (I3.6p - (I3.7I) (i.e., without integrating the Raychaudhuri 

equation). Using 6 = ViU\ we have 

2m 2m 

= r = T , 3.17 

t^r ^/t[(E2 - 4)r2 + 4mr] 

In view of the solutions of 6 obtained earlier in this section by integrating the Raychaudhuri 

equation, the expression in (I3.17P deserves attention. This expression for 6 shows that a 

caustic forms at a turning point of the geodesic motion (i.e., where f = 0). For r O"*" 

(r 0^), we have focusing (defocusing). In order to obtain the explicit A dependence of 9, 

one has to substitute r(A) from the solution of the geodesic equation for r in (I3.17p . It may be 

easily checked that the solutions thus obtained from (13.170 are special cases of the previous 

solutions. It is important to note that, unlike the solutions of 6 obtained by integrating the 

Raychaudhuri equation, in the expression (I3.17p . there is no way of specifying any initial 

condition on 6. Thus, using (13.170 . one cannot study the effect of initial conditions on the 

evolution of a geodesic congruence. This is a subtle issue which will be discussed further in 

the following section. 



D. Analysis of geodesic focusing 



From the exact solutions of the expansion scalar obtained above by integrating the Ray- 
chauduri equations, one can determine the occurrence of finite time singularity (i.e., caustic 
formation). For Case (A), it may be deduced from (13.91) that the expansion scalar 6 ±00 
as A — s> 7r/2. Thus, we may have focusing or defocusing of geodesic congruences depending 
on the initial conditions. One may therefore calculate a critical initial value of the expan- 
sion scalar by exploiting the indefiniteness condition on 6 which leads to Di = —tt/2 ( for 
Ao < ^/2). Now from (13.101) . one can calculate the critical initial value of the expansion 
scalar as. 



From this analysis and using Aq = 1, we have 6^ ~ —1.95. Thus, for 6q < 6q, we have the 
congruence focusing, i.e., a finite time singularity occurs. 

In Case (B), it may easily be concluded that caustic in the geodesic congruence forms 
for any initial condition 6q < — I/Aq. In this case, we have geodesic focusing whenever the 
initial condition satisfies this condition. 

For Case (C), it may be observed from the solution (13.141) that defocusing is not possible. 
However, focusing can occur for an appropriate choice of initial conditions which can be 
obtained by choosing 60 < 6q where 6q = —2/ sinh2Ao is the critical value of the initial 
expansion scalar. For ^0 > finite time singularity cannot occur. 



It is well-known from the work of Tipler [23| that if RimU^'u™' > (timelike convergence 
condition) then focusing (and conjugate points) arise in the congruence within a finite value 
of the affine parameter. From our above analysis, it may thus seem counter-intuitive that 
Rimvf'u™' = — -R/2 < leads to focusing of timelike geodesic congruences in two dimensional 
spacetimes. This, however, is not in conflict with the results of Tipler. In situations, 
such as those presented above, the timelike convergence condition is clearly violated. But, 
with appropriate initial conditions on the expansion (as shown above), one may still have 
a focusing of geodesic congruences. Therefore, we may say that initial conditions have a 
crucial role to play in focusing. 

The Case (C) (i.e., with E"^ > 4) brings out a subtle and interesting difference between the 
solutions of 6 in (I3.14p and (I3.17p . As is clear from the expression of the effective potential in 
(13. 7p . all outward trajectories (outside the horizon) escape out to infinity without any turning 




(3.18) 
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point. Hence, one may conclude from (13.170 that there are no caustics in such a scenario. 
However, from the solution fl3.14p discussed above, we do have focusing depending on the 
initial condition 6q. This difference can be reconciled with if we realize that the expression 
( I3.17P actually yields the expansion scalar field corresponding to the (static) velocity vector 
field = {t, r). On the other hand, the expression (13.141) tells us about the expansion 
history of a congruence, which may have been started with an arbitrary initial expansion, 
as observed in the local frame of a freely falling observer. This is also the approach adopted 
while proving the well-known focusing theorem 



IV. KINEMATICS OF DEFORMATIONS IN 3 + 1 DIMENSIONS 
A. The evolution equations 

In four dimensions, for a congruence of time-like geodesies, the transverse metric on a 
space-like hypersurface can be expressed as, 

K(5 = Qap + u^up, (a,/3 = 0, 1,2,3), (4.1) 

where is the time-like vector field tangent to the geodesic at each point satisfying u^u" = 
— 1, and the four dimensional metric gap is defined through the line elements (12.21) and (12.31) . 
The transverse metric satisfies u'^hap = 0, i.e., h^is is orthogonal to u"'. This transverse 
space-like hypersurface represents the local rest frame of a freely falling observer in the 
given spacetime. The point of interest in this investigation is to determine the deformations 
in this local rest frame as perceived by the observer. The evolution of space-like deformations 
on this transverse hypersurface can be quantified using the tensor Bap, which can now be 
decomposed as follows, 

Bal3 = -(^ ha/B + CTa/B + ^^a(3, (4.2) 

where 9 = B°^ is the expansion scalar, while aap = B(^ap) — Ohap/S and = B[ai3] 
are the shear and rotation tensors. The brackets ( ) and [ ] denote symmetrisation and 
antisymmetrisation, respectively. The shear and rotation tensors also satisfy /i"^ cTq^ = 
and h""^ uja(3 = 0, as can be easily checked. The evolution equation for Bap takes the form. 

Bap + Ba^B'< = -Ra^ps u^'u'. (4.3) 
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Using this equation, one can now obtain the evolution equations for the ESR variables which 
are discussed below. 



1. Raychaudhuri equations 

Following well-known methods, the Raychaudhuri equations for the expansion scalar, 
and the shear and rotation tensors can be obtained as, 

+ \9^ + ia^~uj^) + Ra/3u''u'^ = 0, (4.4) 

2 1 1 ~ 

c^afs + ^0 ao,p + cr„^cr'^^ + u^a^u)^/^ ~ 3 ~ ^"^^ ^ ^""""^^ '^^ " 2 ^""^ ^ '^^'^^ 

2 

+ ^0 + CtJ U^p + CuJ 0"7/3 = 0, (4.6) 

where cr^ = aaf3 cr"'^, uj"^ = ujap cu"^, Ra/s = hpsR'^^ — haf3 h^&R^^ /?> and Caj3r,& is the Weyl 
tensor. These equations are first-order, coupled, nonlinear and inhomogeneous differential 
equations. The equation (14. 4p for the expansion is the well-known Riccati equation, and, 
as mentioned before, is of prime importance in the context of the proof of the singularity 
theorems in general relativity 8|, l9| and in establishing the notion of geodesic focusing |6| . 
With the projection metric defined in (14.11) . the Raychaudhuri equations essentially turn out 
to be structurally similar as in the case of three spatial dimensions. It may also be noted 
that there can be some congruences having a vanishing vorticity for which the velocity vector 
field is hypersurface orthogonal, and (14.61) becomes identically zero. 



2. First integrals of geodesic equations on the equatorial section 

One may note that the equations (12. 7p and (12.81) are independent of X(r), Y{r) and their 
derivatives. Without the loss of generality, one can choose = ti/2 which satisfies (12. 7p 
identically. With this choice, we can capture the kinematics of deformations in the r-0 plane. 
One can now integrate (12.81) once to obtain = C/r^ where C is a constant of motion. We 
will hereafter use these considerations. It is also noteworthy that with ip = 7r/2, this four 
dimensional description reduces to a three dimensional one and the Raychaudhuri equations 
for the components of shear and rotation can be calculated in a way similar to our recent 
work (see the reference |iq|). In addition, we must keep in mind that the time-like vector 
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field satisfies the normalisation condition u'^Ui = — 1, which leads to, 

r^[-X{r)P + Y{r)r^ +1]+C^ = 0. (4.7) 

The above constraint (14.71) also represents a first integral of the set of geodesic equations 
(I2.5p -( l2y8l) for a specific choice of the constant of integration. We will now discuss the 
geodesic equations and effective potentials for the cases corresponding to the line elements 
(1^ and (23]) respectively. 



Case I (corresponding to line element (12. 2p ) 



The first integral of equation (12. 5p for this case is calculated as follows, 

E (r + m sinh^ 

t = r . (4-8) 

2 r (r — m) 

Now, using (14. 8 p in the constraint (14. 7p leads to, 

1 4C^ 
= — - [E"^ [r + msinh^ a)^ + 4mr - 4r^ H (m - r) ] = -V^ (r), (4.9) 

where V^{r) is the effective potential for the case of electric black hole. The effect of the 
parameter a on the radial motion for different values of E and C can be visualised directly 
from (14.90 . The orbits (circular, scattering and plunge) appear to be qualitatively similar 



to those in Schwarzschild geometry (see [24|). 



Case II (corresponding to line element (12. 3p ) 

The first integral of equation (12. 5p for the magnetic case is as follows, 

^_E{mr- g2) 
2 m (r — m) 

The constraint (14.70 along with equation (I4.10p then leads to. 



(4.10) 



r^ = -r^ [E^ + — (4mr - + (m - r) ) ] = -Vf (r), (4.11) 

where V^{r) is the effective potential for the magnetic black hole. The structure of orbits 
are qualitatively same as in Case I. 
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(r) 
(r) 
(r) 



FIG. 1: Effective potentials for magnetic (V4 (r), Q = 0.25), Schwarzschild {V^{r)) and electric 
(V2^(r), a = 0.25) black holes with E = 1.95 and C = 2.2. 

The causal structure of the electric and magnetic stringy black hole spacetimes is similar 
to the Schwarzschild geometry [20]. This provides a motivation for comparing these three 
cases. The Schwarzschild metric can be constructed from (12.41) with, X{r) = (1 — m/r) 
and Y{r) = (1 — m/r)^^ where usually m = 2M with M as the mass of the Schwarzschild 
black hole. Later, we will consider m = 1 for numerical computations. The Raychaudhuri 
equations in the Schwarzschild case follow from (I4.4p - fl4.6p with Rapu'^u^ = and Rai3 = 0. 
The geodesic equations are well-known and the first integrals of the t and equations are 
same as those for the stringy black holes in 3 + 1 dimensions (see Section IIV A 2p with 
ip = tt/2. The first integral of the r equation which satisfies the time-like constraint (14 .yp 
on the velocity field leads to the following effective potential. 



4(7^ 

{E'^ — 4)r^ + 4mr H (m — r) 



(4.12) 



Since the Schwarzschild metric reduces to that of the fiat spacetime for r ^ 00, the 
kinematics of deformations far from the singularity is same as for a static fiat spacetime, a 
case which has already been well studied (see [HI]). It may be noticed that the potential (14. 9 p 
(or (14. lip ) reduces to (I4.12p for a = (or Q = 0). The effective potentials are graphically 
presented in Fig. □ (for E = 1.95) and in Fig. [2] (for E = 2.01). 

For the potential in Fig. (H one can have bound or infalling trajectories (depending on 
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FIG. 2: Effective potentials for magnetic {V^^{r), Q = 0.25), Schwarzschild {Vf{r)) and electric 
{Vf{r), a = 0.25) black holes with E = 2.01 and C = 2.2. 

the value of tq), while the potential in Fig. [2] allows infalling and escaping trajectories. 
Such trajectories are observationally important and therefore, in the following, we study the 
kinematics of deformations, numerically, in the above-mentioned backgrounds. 



B. Analysis of deformations in the equatorial section 

In this section, we study the kinematics of deformations restricted to the equatorial 
section of the black hole background. Consider the = 7c/2 section of the spacetime which 
is a 2 + 1 dimensional slice with an induced metric, say, 70,/?. The timelike geodesic motion 
is now confined to the r — (f) plane. The geodesies of the induced metric are also geodesies of 
the full metric, and the effective potentials remain unaltered. We choose the deformations 

to be 2 + 1 dimensional and accordingly we define Ba/3 = V/jWa where = (j;,r,^^. 
Similarly, one can also define hap = 'Jais + UaUp. 

We then consider a freely falling (Fermi) normal frame E'^ (with Ef = 11°") which is 
parallely transported according to u"Va E^ = 0. The kinematics of deformations, restricted 
to the 2-dimensional spacelike hypersurface (representing the local frame of a freely falling 
observer) in this basis, can now be represented by four kinematical quantities, namely, 9, 
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cr+, (Jx and OJ. The tensor 13 in this basis, can be constructed as, 

Bap = i^e + a^)ele'^ + {^6 - a^)eie\ + (a, + oj)ele^^ + (cr^ - oj)e% . (4.13) 

where are co-frame basis satisfying e'^E^ = 5^. The ESR can be extracted from the 
evolution tensor ( 14.13P using the basis vectors as follows, 

^^ = 5,^/i"^ = 5,^7-^ (4.14) 

a+ = ^ (B^f, E^ - E^E^^), (4. 15) 

1 
" 2 
1 
2 

As in the 1 + 1 dimensional example discussed earlier, the first integrals of the geodesic 
equations of the 2 + 1 dimensional line element enable us to find the expansion, shear and 
rotation for a geodesic congruence. Making use of the vector field and the definition of 
Bai3, we can obtain 6 for example as follows, 

1 €-^-ir' + C')i^ 



a.=^ {B^p E^ E^ + B^p E^ E^), (4.16) 
u; = ^ {B^p E^ E^ - E^). (4.17) 



rXY 



^ E"^ X 2 



4r r 2r 



4r r ^ ~ ' 2r^ 



(4.18) 



where X{r) and Y{r) are the metric functions defined earlier. Similar general expressions 
for a+, o"x and can also be obtained. It may be noted that the expressions thus obtained 
are all functions of r. Obtaining r(A) by solving the geodesic equations, one can then find 
0{X), crjj(A) and uJij{X). From all these expressions, it is easy to state that divergences in the 
ESR appear at the turning points (i.e. where f = 0). Similarly as in 1 + 1 dimensions, these 
solutions represent only special solutions with special initial conditions. They do not reveal 
the effect of initial conditions on the evolution of the ESR variables. A more general class 
of solutions corresponding to arbitrary initial conditions on the ESR variables are obtained 
by integrating the full set of Raychaudhuri equations together with the geodesic equations. 

In order to understand the caustic formation/focusing behaviour in more detail, let us 
redefine the expansion scalar as 6* = 2F/F where the dot indicates the derivative with 



respect to A. One may then use the Fermi normal basis to rewrite (14.41) as the following 
Hill-type equation, 

F + HF = 0, (4.19) 
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where H = al + al - uj'^ + G G = m° 

The notion of focusing is related to F = 0, F < at a finite A. This can be achieved 
under specific conditions on (a) the sign of H ^ and (b) the initial values of the ESR. For a 
complete analysis, we would require to consider all possible initial values (positive, negative 
or zero) for the ESR and H . Here, we restrict ourselves to some special cases and briefly 
comment on the rest. 

(A) \{ H > Q for all A, then conjugate points exist and focusing takes place, as is well 
known. 

(B) When if < for all A, we have focusing only when < < 0. When 9q > 6q, there 
is defocusing. 

(C) If H is sign indefinite over the range of A, then there may be various possibilities 
depending on the initial conditions of the ESR. 

The conclusion (B) has already been illustrated above for the 2-dimensional case in Sec- 
tion [ml 

We now illustrate the above conclusions with our numerical evaluations and corresponding 
plots for the ESR variables. In what follows, we analyse and compare the kinematics of 
deformations in the charged (electric and magnetic) stringy black hole with the Schwarzschild 
black hole. For this, we consider the two potentials shown in Figs. [1] and [2] for the bound and 
escaping trajectories for the reason mentioned before. Corresponding to these two potentials, 
the ESR variables for the two black holes are compared below. In both the cases, the initial 
radius is taken as rg = 8. The initial r is considered as positive for the evaluations for Fig. [31 
while f is taken as negative for the evaluations depicted in Fig. [HI 

In Fig. [3| (corresponding to the the potential shown in Fig. [1]), we observe geodesic focusing 
in all the three backgrounds. In order to understand this behaviour in light of the general 
conclusions mentioned earlier, we have plotted H and G in Fig. 4, for geodesic congruences 
in the electric and magnetic stringy black holes. It is evident that G is negative definite 
along the geodesic flow. However, H is initially negative over a certain range of A but soon 
turns positive and increases monotonically due to dominating contributions from the shear 
components. Thus, one may say that the focusing effect seen here is largely shear-induced. 
Even though the rotation term appears with an opposite sign in if (A), the net value, beyond 
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FIG. 3: Comparison of ESR variables in electric (a = 0.25), magnetic {Q = 0.25) and Schwarzschild 
black holes with E = 1.95, C = 2.2, 6*0 = 0.2, ct+q = 0, (ixo = 0, = 0.1 and ro = 8. 
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(a) 



(b) 



FIG. 5: Variation of time to singularity (A^) versus initial expansion (6q) for different black hole 
metrics (a) without initial rotation of the congruence, and (b) with initial rotation of the congruence 
{ujQ = 0.1). Here, we have considered cr+o = and cjxo = 0. 



a certain A becomes positive and may be responsible for the focusing effect. We may note, 
in passing, that the focusing effect in Schwarzschild spacetime is entirely due the eventual 
dominance of shear over rotation. 

The time of approach to a singularity in the congruence, which we denote by As, can 
also be an interesting quantifier which we can use to characterise geodesic flows in the three 
backgrounds. We have numerically studied the effect of the initial expansion 6q on A^ for 
the different black hole metrics without and with initial rotation ujq of the congruence. The 
results of this study are shown in Fig. [51 It is interesting to note that for an initially 
contracting congruence the singularity occurs more rapidly as compared to an initially ex- 
panding congruence. Furthermore, the time to singularity does not change appreciably for 
large values of initial expansion/contraction. On the other hand, even with a small initial 
rotation {ujq = 0.1), A^ remains almost unchanged over the whole range of variation of the 6q 
considered. It was also found (though the results are not presented) that, with initial shear, 
Ao reduces drastically over the complete range of 6*0 which is also expected qualitatively (see 
[61]). One difference that emerges from this study on the three black hole metrics considered 
is that Af > Af > Af^ as observed in Fig. [5l 

In Fig. [6] (corresponding to the the potential shown in Fig. [2]), there is no focusing. We 
observe from Fig. [7] that H is sign indefinite and bounded over the range of A. As mentioned 
in (C) above, in such situations, definite conclusions are difficult to arrive at because of the 



19 



0.4 
0.3 
0.2 
0.1 




/\ 




/ \ 





/ \ 









% 








\ 

















50 



100 



0.08 
0.06 
0.04 
0.02 




;\ 
1 \ 
1 \ 
1 1 
1 \ 


, S\2 

- - (O' 















50 
% 



100 



0.1 
0.08 
0.06 
0.04 
0.02 














■ - - 03^ 












co^ 







10 



20 30 



40 50 



FIG. 6: Comparison of ESR variables in electric (a = 0.25), magnetic {Q = 0.25) and Schwarzschild 
black holes with E = 2.01, C = 2.2, 6*0 = 0.2, ct+q = 0, a^Q = 0, = 0.1 and ro = 8. 

sensitive dependence of the results on the initial values of the ESR. The following curious 
features may be associated with the no-focusing behaviour observed above. Firstly, it is 
easy to note from the Fig. [7] that J H dX has a negative value, unlike that observed in Fig. 
4. Further, the nature of H (and H) in Figs. 4 and [71 for large A, are distinctly different 
and may also be a cause of the no-focusing effect. 

Additionally, the location (in A) of the extrema (maxima/minima) in the ESR shifts to 
larger A as we move from the electric to Schwarzschild to the magnetic solutions. This is 
evident in Fig. 6. If Ag denotes the location of an extremum, from the figure, we can easily 
say that Af > Af > Af . 

Finally, it may be noted from the plots in Fig. 6 corresponding to u that the rotation 
of the congruence is largely similar irrespective of the metric. This is expected since u = 
Uo exp(— J 6dt), and the variation of 6 is almost similar in the three backgrounds. In 
contrast, from this expression of uj, we note that the divergence (to — oo) of 6 is reflected in 
the divergence (to +oo) of a;, as shown in Fig. O 
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FIG. 7: Variation of G and H for (a) electric and, (b) magnetic black holes with the initial 
conditions as for Fig. [H 

V. SUMMARY AND CONCLUSIONS 

In this article, we have investigated the kinematics of timelike geodesic flows in two and 
four dimensional spacetimes representing stringy black holes. We now briefly summarise the 
work done, the conclusions drawn from it and also mention possibilities on future work. 

(i) The exact solutions of the expansion scalar for the different cases in 2D have been 
calculated by solving the corresponding Raychaudhuri equation for the expansion. 
The occurrence of a finite time singularity (i.e., caustic formation/geodesic focusing) 
in each case is then discussed with particular reference to the relation between initial 
conditions and the behaviour of the expansion. 

(ii) The geodesic equations and the Raychaudhuri equations for the ESR are written out 
and solved numerically for timelike geodesic congruences in two different stringy black 
hole spacetimes in four dimensions. 

(iii) We have drawn some general conclusions and made some observations on geodesic 
focusing which we believe are new. In particular, we have demonstrated how differ- 
ent initial conditions on the ESR can affect the occurence of geodesic focusing. We 
also show how focusing can be affected by the variation of H{X). Even in situations 
where the timelike convergence condition is violated, domination of shear can still lead 
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to focusing. Further, we have introduced a new quantity - the time of approach to 
singularity, which may be used to distinguish between geodesic flows in different back- 
grounds. Though not presented in this article, we have observed that, in the presence 
of initial shear (i.e., a+o 7^ and/or a^o 7^ 0), the time of approach to singularity (A^) 
is significantly reduced. 

(iv) In the scenario depicted in Fig. 6, where there is no focusing, we make an attempt 
towards understanding why this happens by analysing the behaviour. We observe here 
that the locations of the extrema show a systematic shift as we move from geodesic 
flows in the electric to the Schwarzschild and then to the magnetic black holes. 

(v) On the whole, in some sense, the stringy nature of the black hole geometry does seem 
to manifest itself in the nature of evolution of the ESR. 

An interesting issue that is still left unanswered in this work in the role of duality of 
the electric and magnetic black hole metrics on the kinematics. The question that might 
be asked is whether the kinematics in these two spacetimes are also dual of one another in 
some sense. It may be tempting to approach this issue by searching relations between the 
parameters a, Q, m, E and C which leaves the kinematics invariant. 

The metrics of the stringy black holes have coordinate singularities at specific values 
of r and hence cannot be extended beyond. Thus, for a more complete description of the 
kinematics of flows, it will be interesting to study geodesic flows using a different, extendable 
coordinate system (viz. the maximally extended Kruskal coordinate system). Besides this, 
our work has so far focused entirely on timelike geodesies in static spacetimes. It would 
be worth studying the nature of null congruences in a similar fashion. A logical next step 
would be to consider geodesic flows in stationary metrics (such as the Kerr black hole and 
its generahsations) . 

Finally, the essential goal behind this work has been to demonstrate that the Ray- 
chaudhuri equations and the geodesic equations can be solved simultaneously to give us a 
complete picture of geodesies and geodesic flows in any given spacetime. Thus, we have a 
viable approach for studying the kinematics of geodesic congruences for any given metric 
which can help us distinguish between spacetimes through the behaviour of trajectories and 
families of trajectories. In the long run, it may be possible to make use of these results in 
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arriving at distinct observable effects in specific gravitational fields. 
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